Mesoscopic fluctuations of Coulomb drag between quasi-ballistic ID— wires 
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Quasiballistic ID quantum wires are known to have a conductance of the order of 2e^//i, with 
small sample-to-sample fluctuations. We present a study of the transconductance G12 of two 
Coulomb-coupled quasiballistic wires, i.e., we consider the Coulomb drag geometry. We show that 
the fluctuations in G12 differ dramatically from those of the diagonal conductance Gu: the fluc- 
tuations are large, and can even exceed the mean value, thus implying a possible reversal of the 
induced drag current. We report extensive numerical simulations elucidating the fluctuations, both 
for correlated and uncorrelated disorder. We also present analytic arguments, which fully account 
for the trends observed numerically. 
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I. INTRODUCTION 

Recent advances in nanotechnology have made trans- 
port studies of quantum wires an active field of study. 
Several fabrication routes have become available: gated 
tvsfo-dimensional electron gases, cleaved layer overgrowth 
techniques, grooved high-index surfaces, or even nan- 
otubes. It is quite feasible that in future's nanoelec- 
tronic components several quantum wires are very closely 
spaced (spacing d is of the order of the inverse screening 
wave- vector), and hence the stray fields due to moving 
charges in one quantum wire affect the motion of charges 
in the neighboring wires. Analogous effects in coupled 
quantum wells have been intensively studied for the last 
ten years (for a review, see Ref. |l|), and have become 
known as Coulomb drag. While only very few experimen- 
tal studies of coupled quantum wires have been reported 
so far (however, see Refs. we believe that more ex- 
periments will become available in the near future. 

Coulomb drag of mesoscopic structures has been ad- 
dressed theoretically in the case of lE),|Systems both 
within the Boltzmann equation approachci El and jfoLtiJjut- 
tinger liquids with strong interwire interactions .Lrtll For 
the latter case, the interesting possibility of a regime 
with alm^t. identical currents in the two wires has been 
predictedOij and also interesting effects for drag between 
carbon nanotubes have been found.BEj 

Disorder will inevitably be present in all real samples. 
The study of ffuctuation phenomena in mesoscopic sys- 
tems is a mature field, and has tremendously increased 
our understanding of the basic physics governing elec- 
tronic transport in solids. The observatioui-and expla- 
nation of universal conductance fluctuationx3 is one of 



the central achievements in this field. It is then natu- 
ral to ask: are similar phenomena present in the case of 
Coulomb coupled systems? In other words, what are the 
fiuctuation properties of the transconductance? 
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FIG. 1. (a) Geometry of the two coupled ID-wires of 
length L and separation d. (b) Lattice model of the two 
ID-wires. Here, • denote the lattice points of the wires (where 
Coulomb interaction U12 is effective) and o denote those be- 
longing to the ideal leads. 

Quite recently the study of fluctuations of the jGpulomb 
drag was initiated by Narozhny and Aleiner.llj These 
ideas were used by us to study various properties of 
Coulonjib. drag in systems smaller than the Thouless 
energy.oEj They provided answers to some aspects to 
the question phrased above. In particular, in Ref. |l3| it 
was shown that the fluctuations will be pronounced for 
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temperatures smaller than the Thouless energj/j— In con- 
trast to the universal conductance fluctuationsc£l of the 
Landauer conductance Gu the fluctuations of 6*21 are 
however non-universal and depend on temperature and 
various system parameters. 

In this work we focus on coupled quasi-ballistic ID 
wires (see Fig. |l|) and show that even weak disorder can 
change the ballistic properties drastically, and give rise 
to new interesting phenomena. The paper also gives a 
more detailed derivation of some technical results used 
in our recent Letter. til 

We consider temperatures smaller than hvp/L so that 
intra-wire interactions are effectively given by the Fermi 
liquid reservoirs and where the screening can be consid- 
ered to be static. In this regime it is well known from 
calculation of drag between clean Luttinger liquids tl 
the intra-wire interactions have no pronounced effect! 
For the inter-wire interaction we consider the physi- 
cal relevant regime with longe-range interaction where 

C/l2(g = 0)»C/i2(<Z = 2fc;^). 

Starting from the Kubo formula we derive an expres- 
sion f(Mi.|tte transconductance 021 similar to previous 
resultal3'tZ -but relaxing the assumption of translational 
invariance.Eil We formulate the problem with the help 
of the spectral function and subsequently map it onto a 
tight-binding-like model suitable for a computer imple- 
mentation. Our formalism allows us to treat wires with 
electrons propagating in arbitrary potentials. The fluc- 
tuations in drag between disordered wires is studied by 
numerical ensemble averaging and also analytically with 
the disorder included perturbatively. We consider two 
situations: uncorrelated and correlated disorder. For the 
situation where the disorder potentials of the two-azi 
are mutually uncorrelated, as is usually assumed,E3iii 
we find large fluctuations of the order of the mean value 
and thereby also a possible sign reversal of the drag cur- 
rent. If the two disorder potentials are identical, i.e. 
mutually correlated,t3 we get large fluctuations, and an 
enhanced mean value compared to uncorrelated disorder. 
We also give some results for the intermediate case. The 
case of fully correlated disorder is in qualitative agree- 
ment with recent predictions for 2D systemsll3 and for 
the fluctuations we predict an enhancement by a factor 
of which is confirmed numerically. 

The obtained distributions are quite robust in the sense 
that systems with different disorder strength and/or 
length and separation can be rescaled to fall on the same 
curve. However, the distributions also depend on the 
range of the interaction U12 and in that sense they are 
non-universal. 

The paper is organized as follows: in Sec. |l| we in- 
troduce the formalism and in Sec. [II we derive analytic 
results for quasi-ballistic wires. Sec. [V contains our nu- 



merical results. Finally, in Sec. ^ discussions and con- 
clusions are given. Certain technical details can be found 
in appendices. 
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FIG. 2. Diagrammatic representation of the current-cur- 
rent correlation function Ai)7i2L''i2A2, see Eq. (eI). 



II. FORMALISM 

A. General formulation 

Let us consider two ID-wires of length L (shorter than 
the phase-coherence length parallel to each other 
with a separation c?, see Fig. y. Our aim is to calcu- 
late the dc [Vl — > 0) linear-response transconductance 
6*21 = 012/ dVi. In the ID situation the transcon- 
ductance is equal to the transconductivity which can 
be calculated from Kubo formalisnxJ which expresses 
G21 in tetms of the retarded current-current correlation 
functiontZl 

2 j-co 

G21 = lini — / d(t - t') 

X exp [in{t - t')] { [J2ix2,t), Jiixi,t')]). (1) 

Here, Ji is the particle current operator of wire i = 1,2. 
The brackets indicate the quantum mechanical statistical 
average. Due to current conservation arf_and X2 can be 
chosen at any position along the wires.E3Eil Using Mat- 
subara formalism we follow Ref. ^ and calculate G21 to 
second order in the interaction U12 between the meso- 
scopic ID-wires. Intrawire interactions are neglected, 
except for disorder potential scattering. The expansion 
of the time-development operator to second order in U12 
gives rise to a current-current correlation function which 
can be expressed as a product of two three-point cor- 
relation functions: {pipiJi)Ui2Ui2{p2P2J2), where pi is 
the particle density operator. The three-point correlation 
functions are evaluated with the aid of Wick's theorem 
and after some lengthy, but in principle straight forward 
calculations, we arrive at 
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6*21 



L pL pL nL 



h Jo Jo Jo Jo 
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dx'i dx'2 dx'(dx2 



hduj 



Ai(^, x[,x'l)Ui2ix[;x'2)Ui2{x'{; 4)A2(a;, x'^, 4) 
2kTsmh^{hLj/2kT) 



X \nF{£ — 'huj) — nF{£)\ — { lo ^ —lo )*, 
Alix^x') = i \giix,x') - {giix',x)y] = (x|a)(a|x')(5(£ - 



(2) 



(3) 
(4) 



Here, Al is the spectral function, Ql the retarded Green 
function, npis) the Fermi-Dirac distribution function, 
and Sq is the eigenvalue of the exact single-particle eigen- 
state la) in the uncoupled wires. We furthermore assume 
spin-degeneracy. The current-current correlation func- 
tion AiJ7i2t/i2A2 is shown schematically in Fig. |2[ 

Eq. (||) generalizes the results of Ref. |l^ to broken 
translational invariance, and App. ^ establishes a cop=. 
nection to the expression employed in our recent letter .l3 



B. Lattice formulation 



as sums 



Writing the four spatial integrals in Eq. 
over the functions on a discrete lattice, see Fig. 0, we get 



hduj 



Tr[L/f2Ai(u;)£/i2A2(c^)] 
2kTsmh^{hLj/2kT) 

where the matrices have elements 



G21 — -r 
h 



n'a), 
J de{AUj^,^ 



xTjAn^ ,{A\\ ^^-{A\\ AA\\^^^ , 

L L c J ji-xi' y £ ) n,n+l L ^ J nn L ^ J n~\~l,n' 

x[np{e — Tiuj) — nF{e)\ — ( u; ^ —uj 
{Al\ , ^ ax AUx na,x' — > n'a). 



(5) 



(6) 



(7) 
(8) 



where n,n' — 1, 2, 3, ... TV label the lattice points and a 
is the lattice constant, see Fig. |^(b). In the matrix no- 
tation ^ i[Ql~ [Ql]'^]- The derivative d^^ in A has 
been accounted for by the method of finite differencesEHl 
and n can be any lattice points 1, 2, 3, . . . N — 2, N— 1 due 
to current conservation. Summing over the first iV — 1 
lattice points and dividing by — 1 we find 



The next step is to calculate the lattice representa- 
tion of the retarded Green functions and of the 
two uncoupled wires. Writing the Laplacian in the 
wire Hamiltonians with the help of finite differences the 
problem of the uncoupled-wires is mapped onto tight- 
binding-like Hamiltonian£3 



H^' = (2< + {V^}n)Sna' " tSn,n' ±1 , t = /2ma\ (11) 

A standard approach based on Dyson's equation then 
gives the retarded Green functions of the uncoupled wires 
as an N X N matrixE3 



(12) 



where the couplings to lead p = L, R are taken care of by 
the retarded self-energy 

{E;(e)}„„, = -iexp(z/c(e)a)<5„,„^<5„^,„- (13) 

with riL = 1 and riR = A^. The wave vector is related to 
the energy through the usual cosine dispersion relation 
with a band-width of At, i.e. e = 2t{l — cos ka). 



C. Low temperature expansion 

We now consider the low temperature limit kT <C 
where we can evaluate the spectral functions in Eq. (^ at 
the Fermi level. Performing the energy integration gives 
Ai (X u! and the oj integration in Eq. (0) can now be done: 



hduj ■ 



2fcTsinh2(?iw/2fcT) 
and we get 021 oc T^. Eq. (0) now simplifies to 



(14) 



X [^^(e — huj) — nF{e)] — ( w ^ —cu )*, (9) 

where (E" denotes an element-by-element multiplication, 
{X (g) Y}nm = XnmYnm and the matrix A has elements 



A. 



±5^ 



n,n'±l 



N-l 



(10) 



e 2 I 



G21 = J {kTY -Tr [?7i2 A'fi C/12 M2] , (15) 



where 



M,=Re{{^^^}^0[^^^AA^J}. 



(16) 



Eq. ( p^ forms the basis for all subsequent numerical 
work. 
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The Landauer conductance Gu = dli/dVi of|-the indi- 
vidual wires can be expressed in a similar forn£3 



Gu = ^Tr[Tl{si.)gi^rl{eF){gi^y], (17) 
where the leads are described by 

^lie) = ^[me)-{ms)}^]. (18) 

The matrix formulation is readily implemented on a 
computer and the accuracy can be increased -simply by 
having more lattice points N (for a given L)E3 The re- 
tarded Green function can be obtained either by a di- 



rect inversion as indicated in Eq. (12) or by a recursive 
method, 



see e.g. Ref. g3|. 
Alternatively one can also view Eqs. (^5|) and (^^ as 
formulas for a tight-binding system where t is a 'hopping 
matrix element' between different orbitals and the local 
potential Vn + 2t is the energy of the orbital localized at 
site n. Finally we note that the form of Eq. (|l^) is valid 
for any number of transverse channels (for 2D and 3D sys- 
tems A in Eq. (10) should be modified to also contain an 



'integration' over the transverse direction in order to get 
the total current) and thus provides a versatile starting 
point for investigations of drag in different geometries. 




FIG. 3. Example of a first order connected diagrams in the 
diagrammatic expansion for the fluctuations giving rise to a 
{[SG2i]^) oc (X l/kpi dependence. Due to momentum 
conservation the impurity line must carry g = and the four 
interactionlines 2kF. 



A. Fluctuations of the order of V 

An example of a lowest order connected fluctuation 
diagram is shown in Fig. ^ Due to momentum conser- 
vation the impurity line must carry q ~ Q (corresponding 
to forward-scattering) and the four interactionlines 2k p 
so that we estimate the magnitude as follows 



D. Anderson model 

For the study of disordered-syires we use the Anderson 
model with diagonal disordered where the site potentials 
of a given wire are statistically independent with each 
site energy taken from a uniform distribution of width 
W and zero mean (see also Ref. |2^) . The back-scattering 
mean free path £ can be related to the disorder strength 
W and in the Born approximation we get (see App. IbI) 



{[5G21?) o,V\Q)Ut2{2kF). (20) 

where V{q) and Ui2{q) are the Fourier transforms of the 
disorder potential and the Coulomb interaction potential, 
respectively. Though diagrams of order are paramet- 
rically smaller, we shall see that they give the dominant 
contribution in the case of long-range Coulomb interac- 
tion. 



B. Fluctuations of the order of V"^ 



12a{4teF - sl)/W^. 



(19) 



We will examine several different ways of realizing the 
disorder in the two wires. These include a situation where 
both wires are disordered, but with Vi and V2 mutually 
uncorrelated (as in Ref. 14). We also present result for 



disordered wires where the disorder is fully correlated, 
i.e. Vi = V2 (as suggested in Ref. p^. 



III. ANALYTICAL RESULTS 



In this section we consider quasi-ballistic wires. For 
weak disorder it is possible to make a perturbation expan- 
sion for the fluctuationstj SG21 = G21 — {G21) (brackets 
indicate ensemble averaging). 



Fig. m shows examples of connected diagrams of the 
order of V^; now the impurity lines can carry both q = 
and q = 2kp corresponding to back-scattering. For back- 
scattering in the disorder channel we in the case of un- 
correlated (uc) disorder get the estimates 



([5G2i]2)^^ « V*i2kF)U*^{2kp), 



where to lowest order V^{2kp) cx l/kp£ so that 

{[SG2iry^' 0,1/kFi. 



(21a) 
(21b) 
(21c) 
(21d) 



(22) 
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FIG. 4. Examples of second order connected diagrams in 
the diagrammatic expansion for the fluctuations giving rise 
to a ([5G2i]^) (X (X dependence. Due to mo- 

mentum conservation diagrams with only one impurity line 
do not contribute to the drag fluctuations. The diagrams (a) 
and (c) are relevant for both correlated and uncorrelated dis- 
order whereas the diagrams in panels (b) and (d) are relevant 
for correlated disorder only. 

For a long-range Coulomb potential C/i2(0) ^ 
Ui2{2kp) the dominant contribution comes from diagram 
c which also dominates parametrically over the contri- 
bution to 2nd order, see Eq. (EQ). Screening mainly 
affects the small-g limit of Ui2(q)- The mean value 
{G21) oc Uf2{q = 2fc_F) is thus only weakly affected by 
screening whereas the ^'^-contributions to the fluctua- 
tions will be strongly suppressed due to the presence of 
also Ul2iQ = 0), see Eqs. (^). For sufficiently short- 
range interaction the V^-contributions will eventually 
give the dominant contribution, see Eq. ([2^), but then the 
overall magnitude of the fluctuations will also be small 
compared to the mean value. 

For correlated disorder (c) all diagrams a, b, c and 
d contribute (equally) whereas for uncorrelated disorder 
only the diagrams a and c are relevant. More generally, 
for each topologically different diagram contributing in 
the case of uncorrelated disorder there are two similar 
diagrams contributing equally in case of correlated disor- 
der. Of course there are also other possible diagrams in 
case of correlated disorder but to lowest order in l/kp^ 
this means that 



1/2 



'(a) + (5) + (c) + (d) 



(a) + (c) 



/2 X (a) + 2 X (c) 



(23) 



V («) + (^) 
where (a), (6), (c), and (d) refer symbolically to the 



diagrams in Fig. ^ averaged over disorder. This simple 
argument suggests that the fluctuations in the case of 
correlated disorder will be enhanced by a factor of \/2 
compared to the case of uncorrelated disorder. 

We now turn into a quantitative evaluation of the fluc- 
tuations. In App. ^ we show how Ai(x,x') can be ex- 
pressed in terms of scattering states V'lt {^) ■ We introduce 
a short-hand notation 



G21 - C / C/i2(a:,2/)C/i2(a:'2/')Ai(a:,x')A2(2/',2/), (24) 



with 



A(x,x') =Im{/5+_( 



0}, 



(25) 



where J means an integral over all spatial degrees of free- 
dom (in this case cc, x', y, and y'). Here, p are the par- 
ticle density matrix elements taken in a new basis ^± {x) 
related to the original basis V'± {x) by a unitary transfor- 
mation described in App. The frequency integration 
has been performed to give the T^-dependence and all 
prefactors are included in C, i.e. C oc {e^ /h){kTY . 

Before calculating the prefactor in Eq. (22) in the case 



of uncorrelated disorder we first consider the transcon- 
ductance in the case of ideal ballistic wires. Using the 
free plane- waves ^^±{x) = 'il'±{x) = e^^'^f'^ we get 



A(x, x') — Im{f 



} = sin2fcF(a;' - x). (26) 



The electrons in each wire can be either forward-scattered 
(zero momentum) or back-scattered {2k p) and because of 
momentum conservation both the interaction lines in the 
diagram in Fig. ^ carry momentum 2k p (zero momen- 
tum from forward scattering does not contribute to drag 
in this situation) so that 



G2i{^)=C-UU2kF), 



(27) 



where 

Ui2{q) 



L pL 



dxidx2e*«("i""^'t/i2(a;i,a;2), (28) 



In what follows we often express the results for disor- 
dered systems normalized with respect to the ballistic 
transconductance Gi2(cxd). 

The situation is very different even for weak disorder, 
and the broken translational invariance allows for trans- 
ferred momentum others than 2kp. For the diagrams 
in Fig. ^ the impurity lines carry momentum 2kp cor- 
responding to back-scattering within the wires (forward- 
scattering will not contribute to the fluctuations). How- 
ever, the difference between diagram (a) and (c) is that in 
(a) the interaction lines in each sub-diagram must carry 
the same momentum whereas in (c) one of them can 
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carry e.g. zero momentum while the other carries 2kp. 
Now, since J7i2(0) » Ui2{2kF) for long-ranged interac- 
tion (typically, Uuiq) decays exponentially) the diagram 
(c) gives the major contribution to the fluctuations. This 
can be tested numerically by noting that the diagram (a) 
is relevant to the case where wire 1 is disordered and 
wire 2 is either disordered or ballistic whereas the dia- 
gram (c) is only relevant to the case of both wires being 
disordered. Indeed, by numerically calculating the fluc- 
tuations for a system where one of the wires is ballistic 
and the other disordered we have found a very dramatic 
reduction (by more than an order of magnitude) of the 
fluctuations compared to the case of both wires being 
disordered. 

Considering wires with mutually uncorrelated disorder 
and writing A{x,x') = A{x,x') + SA{x,x') (where the 
second term is the correction to the ballistic limit) the 
contribution from diagram (c) is 

{[5G2if)^^^C^ j Ui2{x,y)Uu{x',y')Uu{x,y) 

xUi2ix' ,y'){SAi{x,x')SAi{x,x')) 
(<5A2(y',y)<5A2(y',y)), (29) 

and the remaining problem is to calculate the correla- 
tion (^SAi{x,x')SAi{x,x')') to lowest order (second or- 
der) in the disorder st reng th. For quasi-ballistic wires 
(T = 1 - 7?. ~ 1) Eqs. (Q and (@ give 



A{x,x') ^lm{p+^{x)p-+{x')}, 



(30) 



with Pajsix) — 'ilj*^{x)ijj p{x) in the original basis and for- 
mally we thus have that 



i2kFX I ^—i2kFX 



5A(x,x') =Im{(5p+_(x)e*2/cj.x _^ 



Sp+-{x')}. 



(31) 

Using the Lippmann-Schwinger equationll we have to 
lowest order in the disorder strength that 



^±(x)~e±^'=-^+ / dx^o(a;,x)y(x)e±^'=-^ 
Jo 



(32) 



where Goix,x') = {iTiVF)~^e-^^''^^ is the unperturbed 

retarded Green function. This means that p^ (x) — 

p'} (x) -|- uix) -\- h'{x) where 



instead of I7i2(2A:f). This is not the case for v and we 
can therefore omit it. For the averages over disorder we 
use {V{x)V{x')) = Wq6{x — x') ^'^'^ neglecting terms 
oscillating with Akp we get (u(a;)u(i)) w and 

{u{x)u*{x)) = L (i^)' (1 - 2\x - x\/L). 

The prefactor can be related to the reflection coefficient 
of the wire. To see this we consider the Dyson equation 
to second order in the disorder 

g[x,x')~go{x,x')+ [ dxgoix,x)v{x)Go{x,x') 

Jo 

dx / dx' go{x,x)V{x)Qa{x,K')V{^)Qo{x' ,x'), 
Jo 

(35) 

and from the Fisher-Lee relation we notice that TZ = 
1 - r = 1 - {hvFy\g{L,0)\^. It then follows that to 
second order in the disorder (7^) ~ L (WoMwf)^ = L/£ 
(the last equality follows from Appendix P) so that 



{u{x)u*{x)) = (7^)(l - 2\x ~ x\/L). 
This means that 



(36) 



{SA{x,x')dA{x,x')) 



(n) 



{f{x' - x',x- x) 



+f{x — x,x' — x') 
—f{x' — x,x — x') 
-f{x - x',x' - x)}, 



(37) 



where f{x'~x',x — x) = {l — 2\x' — x'\/ L) cos2kF{x — x). 

Out of the 16 only the 4 terms of the form cos(a; — 
x) cos(j/ — y) are finite for long wires kpL ^ 1 and we 
thus get 

([^G2l]2)„^~C2(7^l)(7^2) I Uu{x,y)Ui2{x',y') 

xUi2{x,y)Ui2{x',y') 

xf{x ~'X,x' ~ x')f{y -y,y' - y') 

= An^){n2)uU2kF)iJU^). (38) 



This corresponds to the estimate in Eq. (21c). Finally, 
the mean value is close the value for ballistic wires and 
normalizing by Eq. (p7|) we get 



{x) = {invFr^ r dxsigi^{x-x)V{x)p\M). (33) ([^G2i]^)y/ ^ [2(7^l)(7^2)£/f2(2fcf )£/i^ 



,1/2 



v{x) = {ihvF)-^ dxsign{x-x)V{x)p\-{x). (34) 



Here, p\_{x) = exp{—2ikFx). Thus we see that u gives a 
non-oscillating correction to p^_{x) which is the reason 
why the Fourier component J7i2(0) enters the end-result 



G2i(oo) C//2(2fcF) 
where 

pL flj 

Ui2{0)= / / / / dxidx2dx[dx2 
Jo Jo Jo Jo 



(39) 



xU^2{xuX2)UMx\,x',){l - (l - teU), 
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and (Ti-i) = L/£. Eq. ( |39| ) predicts the \/l dependence 
found in the quahtative discussion leading to Eq. (22). 
It follows that the relative magnitude of the fluctuations 
is of the order (TZ')Ui2{Q) /Ui2{2kF) and even though 
(TZ) <C 1 for quasi-ballistic wires the effect of long-ranged 
interaction C/i2(0)/C/i2(2fci?) » 1 can give rise to relative 
fluctuations of order unity, i.e. fluctuations comparable 
to the mean value. As we shall see in the next section, 
Eq. (^9|) is in excellent agreement with our numerical re- 
sults. 



IV. NUMERICAL RESULTS 

In this section we apply Eq. ( [fsl ) to numerically eval- 
uate the statistical properties of ensembles of different 
disorder configurations. 

The numerical implementation was tested by first cal- 
culating the Landauer conductance Ga, Eq. (|l7|), for 
various geometries which can also be solved analytically. 
We have also tested our use of the Anderson model for 
disordered wires. By changing the ratio i/ L from the 
localized regime {i/L <C 1) to the de- localized regime 
{(■/L ^ 1) we have obtained distributions, mean val- 
ues, and fluctuations for the Landauer conductance Ga 
which are iii agreement with the analytical results of 
Abrikosov.E3 We have also compared our lattice imple- 
mentation, Eq. (^5|) to analytical results, see App. for 
drag between two ballistic wires. By making the grid 
sufhciently fine Eq. (|l^ is capable of reproducing the 
curves in Fig. |^ which are based on Eq. (C3) with the 
integrals evaluated numerically. Another test is the pre- 
dicted peakiMjOf G21 at the onset of modes in either of 
the two wirealj and for systems with resonance transmis- 
sion (peaks in Go); we find this to be born out by our 
lattice implementation. 

For the study of drag between disordered wires we 
consider quarter-filled bands [ep = t) and wires with 
TV = 100 lattice points so that fc^i = (tt/S) x 100. The 
separation is kpd = \ and for simplicity we assume un- 
screened Coulomb interaction. 



•1 01234567 



(40) 



For a discussion on how to include an effective screened 
interaction {e.g. within the random-phase approxima- 
tion) see e.g. Ref. |l^ and references therein. 

For the case ^ — 36L (this corresponds ioW — e_F/10) 
the disorder has as expectednZl almost no effect on the 
Landauer conductance, i.e. (Gii) ^ 2e^//i with vanish- 
ing fluctuations. However, for the transconductance the 
situation is very different. Panel (a) of Fig. ^ shows a typ- 
ical histogram of G2i(^)/G2i(oo) (where 6*21(00) is the 
result in the ballistic regime, see App. |c|) for £ = 36L in 
the case of uncorrelated disorder. Depending on the dis- 
order configuration G2i{() can be either higher or lower 




"n~r 

•1 1 2 3 4 5 

G2l{i)/G2l{00) 



6 7 



FIG. 5. Histograms for G2i{£) normalized by the ballistic 
result 6*21(00). The histograms are based on > 10* disor- 
der configurations and for all three histograms the mean free 
path is £ — 36Z/. Panel (a) is for the situation of mutually 
un-correlated disorder, panel (b) for partly correlated disor- 
der (Wc = Wuc), and panel (c) is for mutually fully correlated 
disorder. 



than in the ballistic regime. The enhancement occur- 
ring for certain disorder configurations can be understood 
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physically as follows. The lack of translational invari- 
ance allows forward scattering (transferred momentum 
q ~ 0), which normally has little effect, to cause transi- 
tions between scattering states with opposite directions, 
thus contributing to the drag. We emphasize that the his- 
togram peaks close to the ballistic value and not at zero 
drag, i.e. the mean drag is finite and positive. The vari- 
ance is of the same order as the mean value so that sign 
reversal for some disorder realizations is possible. The 
later is represented by the negative tail in the histogram. 
The sign of the drag is thus arbitrary in the sense that 
both positive and negative drag can be observed. 

For the distribution of G21 we find that Puc(a^) oc 
exp[— — x)/x\°'] with a ~ 1.4 fits surprisingly well 
to the data. However, this observation does not exclude 
the uossibility that other functions may be fitted equally 
well.L^I In fact we have performed these fits to histograms 
for kp£ ranging from 10'^ to 10^ with kpL in the range 100 
to 300 and find indications that a ~ 1.4 — 1.5 which by 
rescaling of G21 makes it possible to let all histograms 
fall onto the same curve. However, the distribution is 
non-universal in the sense that it depends on the range 
of the interaction U12. 



• correlated 




H 1 1 1 — I ■ I I I 1 1 1 1 — I I I I I 1 1 1 1 — I 

FIG. 6. Plot of fluctuations ([5G2i(^)]^)^''^ as a function of 
the mean free path kp£. In both cases the expected l/£ behav- 
ior, Eq. (p^), is born out by the numerical calculations. The 
full lines are Eq. ( P9[) w ith no free parameters. The expected 
enhancement, Eq. for correlated disorder by a factor of 

\/2 compared to uncorrelated disorder is also conflrmed by 
the numerical calculations. 

For the same system parameters but now with fully 
correlated disorder we get a very different distribution as 
seen in panel (c) of Fig. |[ In this situation the mean 
transconductance is enhanced compared to the uncorre- 
lated case. This confirms at least qualitatively the pre- 
dictions for 2D systems of Ref . |l^. The mean fluctuations 
are also enhanced (we return to the enhancement below) . 
However, since the mean transconductance is enhanced 
by almost a factor of two compared to the ballistic limit 
this also means that there is no disorder configurations 



giving rise to negative drag. For partly correlated dis- 
order, panel (b), the distribution is shifted to a higher 
mean value and the shape is also slightly changed with 
an increase in the magnitude of the fluctuations. 

In Fig. H we show the dependence of the fluctuations 
on the mean free path £ which has the expected l/kp£ 
dependence, see Eq. (p2|). We also notice that correlated 
disorder give rise to slightly larger fluctuations compared 
to uncorrelated disorder. In fact they are exactly en- 
hanced by a factor of a/2 as predicted, see Eq. (|2^). 

V. DISCUSSION AND CONCLUSION 

We have studied drag of Coulomb coupled disordered 
mesoscopic ID-wires and developed a formalism for in- 
vestigating the statistical properties both analytically 
and numerically. In this work we have focused on the 
quasi-ballistic regime with L <^ £. For the ordinary Lan- 
daucr conductance Gu this is experimentally a some- 
what trivial regime since weak disorder has almost no 
effect and {Ga) ~ jh with vanishing fluctuations. 
However, surprisingly this same regime (which can be 
achieved experimentally) offers rich physics when mea- 
suring transconductance between two Coulomb coupled 
wires. In this work we report results on the statistical 
properties of the transconductance G21 including its dis- 
tribution function. Depending on the disorder config- 
uration the transconductance can be either positive or 
negative and even though the mean value is close to the 
ballistic value the fluctuations can be of the same order 
of magnitude. This is fully explained by analytical calcu- 
lations including the weak disorder perturbatively. The 
effect can be explained by a combination of disorder in- 
duced back-scattering in the wires and forward scattering 
between the wires induced by Coulomb interaction. 

We have also studied the recently proposed situationEfl 
where the electrons in the two wires experience a com- 
mon disorder potential. In this case we also flnd that even 
weak disorder gives rise to pronounced fluctuations and 
also a considerable enhancement of the mean transcon- 
ductance. 

The reported distribution functions have been success- 
fully fitted by simple analytical expressions. Though 
there are no analytical predictions for the exact func- 
tional form these initial results suggest that a more de- 
tailed analytical investigation could be rewarding. 

Finally, let us note that while this work has focused 
on L ^ i!, for naulti-mode wires also the diffusive and lo- 
calized regimestaEa with J-h> £ seems to be a promising 
direction for future worko 
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APPENDIX A: FORMULATION IN TERMS OF 
SPECTRAL FUNCTIONS 

We start from Eq. (2) of Ref. With a knowledge of 
the exact eigenstates of the isolated wires the expression 
points out a direct way of calculating the transconduc- 
tance and it is possible to simplify this expression for 
ballistic wires where the eigenstates are known analyti- 
cally, see App. ^. However, in most cases the eigenstates 
are not known and it is useful to make a numerically more 
appropriate formulation in terms of the retarded lattice 
Green function. Using that the current matrix element 
can be written as 



{a\j(x)\f3) 



2mi 



lim dx [(a|a;)(5;|/3) 



(a|i)(a;|/3)] , 
(Al) 



and by introducing the factor 1 = J de6{£ — Sa) into 
Eq. (2) of Ref. |l| we get the result in Eq. (|). In the 
second term we have used that the interchange x' ^ x" 
corresponds to a sign-change along with a complex con- 
jugation. 



APPENDIX C: BALLISTIC REGIME 

For ballistic wires the eigenstates are of the form 
'>Pkix) — L~'^l'^ exp{ikx) with e — h^k'^/2m. The spectral 
function, Eq. (0), is then given by 



A^{x, x) — k{e) cos[fc(£)(a;' — x)]/e. 



(CI) 



We consider the low temperature limit kT <C £_f where 
we can evaluate the spectral functions in Eq. (^) at the 
Fermi level so that 



A,(cj,a;',a;") 



2Trh £f 



sin[2fcF ix' ~x")]. 



(C2) 



The u) integration in Eq. can now be performed and 
using Eq. (|lj) we get 



e2 1 /krV (2m> ^ 



<.L rL i-L rL 



pL, pLi p Li pLi 

^ da;^da;2da:'/da:2 sin [2kp {x'l — a;")] 

Jq Ja Jo Jq 

xUi2{x[;x'2)Ui2ix'l; x'^) sin [2k p (4 - 4' )] ■ (C3) 

In general the integrals have to be calculated numeri- 
cally and only in the limit kpL ^ 1 (the limit also studied 
in Ref. ||) we have the asymptotic result 



APPENDIX B: MEAN FREE PATH IN THE 
ANDERSON MODEL 

For a tight-binding chain with N sites and no disorder 
we have eigenstates ipk{n) — N^^^^ exp{ikna). We cal- 
culate the rate for forward (-I-) and back (— ) scattering 
from Fermi's golden rule 



T±{k) 



i±kk' 



(Bl) 



where V is the disorder potential which is treated as a 
perturbation. For a chain with large N we take the sum 
into an integral and since S{£ — e') = |?iufe|~^ [S(k — k') + 
5{k + k')] we get 



1 



T±{k) NTi^ 



\Vk 



nm 



,i{l^l)k{n-m)a rr -xr 



In the Anderson mode0 the different sites are uncor- 
related and p{Vn) = Q{W/2 — Vn)/W which means that 
{VnVm) = {V^)Snni = W^/12 6nm- The Corresponding 
mean free path is given by i{k) = Vkl{Tg^{k)) (forward 
and back scattering give rise to the same result) and for 
£ — 2t{l — cosfca) we get the result in Eq. ( [l9| ) at the 
Fermi level. The result agrees with Ref. |2^ except for a 
constant shift of the energy by 2t. 



ft 6 i £i? 



kT C^(fcF')^^^^^(2/fci.d) 



(C4) 



where Kq is a modified, Bessel function of the sec- 
ond kind of order zeroHj Here, we have assumed un- 
screened Coulomb interaction, see Eq. (^), and intro- 
duced U{r) — e^/(47reo£rr). Fig. shows a numerical 
evaluation of Eq. ( |C3| ) along with the asymptotic result, 
Eq. 



(B2) — 




FIG. 7. Plot of G21 for two ballistic wires, Eq. (|C3|), as a 
function of the length kpL for different values of the sepa- 
ration kpd. The dashed curves show the asymptotic result, 
Eq. (|c3). 
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APPENDIX D: FORMULATION IN TERMS OF 
SCATTERING STATES 

We start from Eq. (4) in Ref. ^ using scattering states 
at the Fermi level as our basis: {'ip+{x),^p-{x)} where + 
(— ) is for a state incident from the left (right) lead. All 
matrix elements Iaf3 and can then be considered as 
elements of 2 x 2 matrices: / and p. Due to current 
conservation we are free to evaluate the current matrix 
element in the leads (outside the region with elastic scat- 
tering) which gives I — {vf/L)J with 

T 

where is the third Pauli matrix and 




(Dl) 



5* = 



r t 
t r' 



-Re'" 



(D2) 



is the usual unitary scattering matrix in the presence of 
time-reversal symmetry [SS'^ — 1 and S ~ S^). In the 
second equality the scattering probabilities T = \t\^ and 
7^ = 1 - r = |rp = |r'p have been introduced. Eq. (4) 
in Ref. 14] can now be written as 



A(c., X, x') = i,^huj^Tr{j[p{x)-p{x')\ ) , (D3) 

where [A; B] — AB — BA is a commutator with A and B 
being matrices. Next we employ a unitary transforma- 
tion 



V* u* 



(D4) 



with |Mp + = 1 that satisfies UJU^ = VTt^ hy 

choosing \u\^ = i(l + VT), \v\^ = i(l - VT), and 

vu* = i-s/^e*'"^"^^ (in a concrete calculation it can 

be useful to use the freedom to choose the phases as 
V = |„|e'(0-e)/2 and u = lule-'^'^-^)/^)^ ^ jg ^j^g^^ g^gy 

to obtain the very compact result 



A(ijj, X, x') — 4:n^huj^^VTilm\^p-^ {x)p ^(x')}, 

(D5) 

where in the new basis p{x) — Wp{x)U. With the choice 
of relative phase mentioned above we have in particular 



p+-{x) = 



T^[p--{x) - P++{x)] 
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